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ABSTRACT 

XMM-Newton observations of 10 Ultra-Luminous Infrared Galaxies (ULlRGs) from a 
200 ksec mini-survey program are reported. The aim is to investigate in hard X-rays a 
complete ULIRG sample selected from the bright IRAS 60 /im catalogue. All sources 
are detected in X-rays, 5 of which for the first time. These observations confirm that 
ULlRGs are intrinsically faint X-rays sources, their observed X-ray luminosities being 
typically L2-iokeV < lO^^"^-^ erg/s, whereas their bolometric (mostly IR) luminosities 
are LhoZ > 10"*^ erg/s. In all sources we find evidence for thermal emission from 
hot plasma with a rather costant temperature kT ~ 0.7 keV, dominating the X-ray 
spectra below 1 keV, and likely associated with a nuclear or circumnuclear starburst. 
This thermal emission appears uncorrelated with the far-lR luminosity, suggesting 
that, in addition to the ongoing rate of star formation, other parameters may also 
affect it. The soft X-ray emission appears to be extended on a scale of ~30 kpc for 
Mkn 231 and IRAS 19254-7245, possible evidence of galactic superwinds. In these two 
sources, in IRAS 20551-4250 and IRAS 23128-5919 we find evidence for the presence 
of hidden AGNs, while a minor AGN contribution may be suspected also in IRAS 
20100-4156. In particular, we have detected a strong (EW~ 2 keV) Fe-K line at 6.4 
keV in the spectrum of IRAS19254-7245 and a weaker one in Mkn 231, suggestive 
of deeply buried AGNs. For the other sources, the X-ray luminosities and spectral 
shapes are consistent with hot thermal plasma and X-ray binary emissions of mainly 
starburst origin. We find that the 2-10 keV luminosities in these sources, most likely 
due to high-mass X-ray binaries, are correlated with LpiR- both luminosities are good 
indicators of the current global star formation rate in the galaxy. The composite nature 
of ULlRGs is then confirmed, with hints for a predominance of the starburst over the 
AGN phenomenon in these objects even when observed in hard X-rays. 

Key words: galaxies: infrared, galaxies: starburst, galaxies: surveys, galaxies: evo- 
lution, galaxies: active 



1 INTRODUCTION 



* Based on observations obtained with XMM-Newton, an ESA 
science missions with instruments and contributions funded by 
ESA Member States and the USA (NASA), 
f e-mail: franceschinini@pd.astro.it 



Ultra-Luminous InfraRed Galaxies (ULlRGs, sources with 
bolometric luminosity Lm > lO^^L© mostly emitted in the 
IR) have received much attention since their discovery dur- 
ing the IRAS survey follow-up observations (Sanders et al. 
1988; see a review in Sanders & Mirabel 1996). The main rea- 
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son for this interest was that, together with optical quasars, 
these are by far the most luminous objects in the local uni- 
verse. 

The relevance of this class of sources in the cosmologi- 
cal context has been further emphasized by recent findings 
of cosmological surveys at long wavelengths. Observations 
at IR and sub-millimeter wavelengths have discovered that 
luminous and ultra-luminous IR galaxies, which are rare in 
the local universe (Soifer et al. 1987), are detected instead 
in large numbers in deep IR surveys, and are a fundamen- 
tal constituent of the high-rcdshift galaxy population (e.g. 
Small et al. 1997; Genzel & Cesaxsky 2000; Franceschini et 
al. 2001). The number and luminosities of these sources im- 
ply that an important fraction of stars in present-day galax- 
ies, or alternatively of the degenerate baryons contained in 
nuclear supermassive Black Holes, have formed during such 
IR-luminous dust-extinguished evolutionary phases in the 
past. It has been argued that luminous and ultra- luminous 
IR galaxies at high redshifts could trace events of star- 
formation which may be at the origin of massive elliptical 
and SO galaxies (e.g. Franceschini et al. 1994; Lilly et al. 
1999). Support for this view has recently come from high- 
resolution spectroscopy of ULIRG mergers (Genzel et al. 
2001), whose dynamical properties have proven that they 
indeed are "ellipticals in formation". 

This evidence came along with the discovery by COBE 
of a bright diffuse radiation, the Cosmic IR Background 
(CIRB, see Hauser et al. 1998; Lagache et al. 2000), ap- 
pearcntly containing a large fraction (up to ^ 70%) of the 
total extragalactic background energy density from radio to 
X-rays, i.e. a major part of the photon energy released by 
cosmic sources at any redshifts. As discussed by various au- 
thors (Elbaz et al. 2002; Small et al. 2002), there is precise 
evidence that this background radiation is indeed produced 
by sources similar in all respects to luminous and ultra- 
luminous IR galaxies at z > 0.5. These recent facts justify 
the growing interest in local ULIRGs as possible clues to 
their supposed high-z counterparts. 

However, about 15 years after their discovery, the na- 
ture of these sources still remains rather enigmatic. Large 
gas and dust column densities in the galaxy cores, respon- 
sible for the IR emission, prevent a direct observation of 
the primary energy source and the IR spectral shapes are 
highly degenerate with respect to the illuminating spectrum. 
It is widely accepted that both starburst and AGN activity 
may be responsible for the observed luminosities (Genzel et 
al. 1998; Veilleux, Kim, Sanders 1999; Risaliti et al. 2000; 
Bassani et al. 2003, in prep.). If indeed ULIRGs are "form- 
ing spheroids", an AGN/starburst association is naturally 
supported by the evidence that all spheroidal galaxies host 
quasar relics in the form of supermassive BH's. The same 
radial inflow of gas produced by the merger is likely to fuel 
both the starburst and, at some stages, the AGN. 

So far conflicting evidence has been reported about 
the relative contributions of the two energy sources in 
ULIRGs, stellar and gravitational. Important progress has 
been achieved with mid- and far-IR spectroscopy, probing 
the inner optically-thick nuclei. Using diagnostics based on 
coronal line intensities and PAH line-to-continuum ratios, 
Lutz et al. (1996), Genzel et al. (1998) and Rigopoulou et 
al. (1999) have argued that the majority of ULIRGs are 
powered by star-formation. 



On the other end, the presence of energetically impor- 
tant obscured AGNs in ULIRGs has been revealed by op- 
tical and near-IR spectroscopy, often detecting Seyfert-like 
nuclear emission line spectra (Sanders et al. 1988; Veilleux, 
Sanders & Kim, 1997) as well as evidence for completely 
buried AGNs (Soifer et al. 2001; Imanishi, Dudley & Mal- 
oney 2001). 

In principle hard X-ray observations offer an additional 

important tool to investigate the presence of hidden AGNs, 
providing quantitative estimates of their contribution to the 
bolometric luminosity. This diagnostics relies on the funda- 
mentally different spectra and luminosity regimes between 
starbursts and AGNs in hard X-rays. Hard continuum emis- 
sion and prominient the Fc-Kq line(s) are distinguishing fea- 
tures of buried AGNs, which can penetrate large gas column 
densities. Some of the brightest nearby ULIRGs, classified 
on the basis of IR spectroscopy as pure starburst, show spec- 
tral properties typical of obscured AGNs when observed in 
hard X-rays (e.g. NGC 6240, Iwasawa et al. 1999; Vignati et 
al.l999). In practice, so far the utilization of this diagnos- 
tics has been limited by the X-ray faintness of most ULIRGs 
(see previous unsuccessful detection attempts with ROSAT, 
ASCA and BeppoSAX; Risaliti et al. 2000). 

The unique large collective area and hard X-ray re- 
sponse of XMM-JVewton have been used here to deeply sur- 
vey for the first time a representative sample of 10 IRAS- 
selected ULIRGs for which high-quality mid-IR and optical 
spectroscopy data are available. Most of these ULIRGs were 
previously classified as starbursts (Lutz, Veilleux & Genzel 
1999). The purpose of the present observations is not only 
to identify AGN signatures in hard X-rays, but also to sep- 
arately quantify the contributions of the buried AGN from 
the starburst emission to the X-ray spectra, and to compare 
them with independent estimates based on optical and IR 
observations. This paper summarizes the main X-ray prop- 
erties of the sample, five sources of which are detected in 
X-rays for the first time. XMM- Newton observations for two 
sources of the sample of particular interest, the Superanten- 
nae and Mkn 231, are discussed in separate papers (Braito 
et al. 2003a and 2003b). 

Section 2 summarizes the properties of our ULIRG sam- 
ple. Section 3 details the XMM-iYcwton observations and 
data analysis. The resulting hard X-ray properties are re- 
ported in Section 4, while Section 5 discusses plausible ori- 
gins of various X-ray components and compares them with 
data at other wavelengths. Section 6 contains our conclu- 
sions. Ho = 50 Km Mpc~^ has been assumed throught. 



2 THE ULIRG SAMPLE 

The sample of 10 ULIRGs discussed in this paper has been 
selected from the list of IRAS ULIRGs observed with ISO 
by Genzel el al. (1998). The sample is flux-limited at 60/im 
and complete to Seo/j > 5.4 Jy. It includes sources with 
L > 10^^ L© in the total 8-1000 nm band, whose far-IR 
selection makes it unbiased with respect to absorption and 
representative of the most luminous galaxies in the local uni- 
verse. High quality IR/optical spectroscopic data are avail- 
able for all the sample sources (Lutz et al. 1999; Veilleux et 
al. 1999). 

Of the original Genzel et al. sample of 15 ULIRGs, 4 
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Figure 1. DSS2 images (2' X 2') of IRAS 12112+0305, IRAS15250 + 3609, IRAS17208 - 0014 and IRAS19254 - 7245 (from top left to 
bottom right). Countours of the X-ray (0.2—10 keV) emission have been overlayed. The countours displayed correspond to 4cr, ba, 7cr, 
10a, 20(7, 30(T, 50(7 above the background. The countours have been derived using the MOS2 data, since this detector has the best point 
spead function (see Ehele et al. 2001). 



have been observed with yiMM- Newton by different teams, 
while IRAS 23060+0505 has been deeply observed with 
ASCA (Brandt et al. 1997) and BeppoSAX. Source names, 
sky positions (which are also the coordinates of the XMM- 
EPIC observations), IR (8-1000 fim) luminosities, redshifts, 
Eis well as the XMM-EPIC net exposure times and observa- 
tion dates, of the 10 ULIRGs presented here, are reported 
in Table 1. Optical and IR classifications of the ULIRGs 
presented here are summarized in Table 4. 



All 10 objects have been observed with XMM-Newton 
with ~20 ksec exposure each. A fraction of the observations 
had to be excluded from the analysis because of background 
flares (this problem was severe for the source IRAS 23128- 
5919, see Table 1). 
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Figure 2. DSS2 images (2' X 2') of IRAS 20100-4156, IRAS20551 - 4250, IRAS22491 - 1808 and MKN231 (from top left to bottom 
right). Countours of the X-ray (0.2—10 keV) emission have been overlayed. The countours displayed (from the MOS2 data) correspond 
to 4(T, 5cr, 7(T, lOcr, 20(7, 30(T, 50a" above the background. See also caption to Fig.l 



3 XMM-NEWTON EPIC DATA PREPARATION 
AND ANALYSIS 



The yiWNl- Newton observations presented here have been 
performed between March 2001 and November 2002 with 
the EPIC (European Photon Imaging Camera: Striider et al. 
2001 and Turner et al. 2001) cameras operating in full-frame 
mode. Data have been processed using the Science Analysis 
Software (SAS version 5.3), and have been analyzed using 
standard software packages (FTOOLS 5.0, XSPEC 11.0). 



The latest calibration files released by the EPIC team have 
been used. 

Event files produced from the pipeline have been fil- 
tered from high-background time intervals, and only events 
corresponding to pattern 0-12 for MOS and pattern 0-4 for 
PN have been used (see the XMM-A^euiion Users' Handbook, 
Ehle et al. 2001). We have then generated our own response 
matrices (that include the correction for the effective area) 
using the SAS tasks arfgen and rmfgen. The screening pro- 
cess from high background time intervals yielded net expo- 
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Figure 3. DSS2 images and XMM-Newton contours (2' X 2') of IRAS 14348-1447 and IRAS 23128-5919. Countours of the X-ray 
(0.2—10 keV) emission have been overlayed. The countours displayed (from the MOS2 data) correspond to 4cr, 5cr, 7cr, lOcr, 20(t, 30cr, 
50cr above the background. See also caption to Fig.l. 



sures between 15.7 ksec and 21.7 ksec (net exposures for PN 
camera are reported in Table 1), with the exception of IRAS 
23128-5919 (9.4 ksec). 

All the 10 ULIRGs have been detected in the MOSl, 
MOS2 and PN detectors with signal to noise ratios S/N>3. 
Close-up images of the sources are reported in Figs. 1, 2, 3 
as overlays of the (0.2-10 keV) X-ray contours on top of the 
POSS red-band images. For all the ULIRGs the bulk of the 
X-ray emission is positionally concident (within the XMM- 
Newton positional error) with the optical-IR core. By ac- 
counting in detail for the EPIC PSF, we have found evidence 
for extended emission only in Mkn 231 and IRAS19254- 
7245, while the surface brightness distributions for all other 
sources are consistent with being unresolved by EPIC. 

Except for Mkn 231, the most extended source (90% 
counts in ~ 40 arcsec in the soft X-ray band E< 2 keV), the 
source spectra were extracted from circular regions of ~ 20 
arcsec radius. Background spectra have been extracted from 
source-free circular regions (with typically ~ 1' radius) close 
to the target. 

In order to improve the statistics, MOSl and M0S2 
data have been combined. Both MOS and PN spectra were 
then automatically rebinned in order to have a > 20 counts 
in each energy channel. The combined MOS data have been 
fitted simultaneously with the PN data, by keeping the rel- 
ative normalizations free. 



4 RESULTS 

We report in Table 4 the measured fluxes and luminosities 
of the sources in various X-ray bands. Our XMM-Newton 
observations show that ULIRGs are rather faint in X-rays if 
compared with AGNs of similar bolometric (IR) luminosi- 



ties. Their observed X-ray luminosities are close to or lower 
than I/2-iofceV = 10*^ erg/ s, while their IR bolometric out- 
put is always Lbo! > 10*^ erg/s. 

The level of detail for our spectral analysis is diflerent 
for each source depending on the qualities of the XMM spec- 
tra, ranging from fairly detailed for stronger sources to only 
rather coarse spectral fits for the faintest. The two brightest 
sources, IRAS 19254-7245 and Mkn 231, have data of suffi- 
ciently high quality to guarantee investigations of both the 
continuum emission and the Fe-K 6.4 keV line, detected only 
in these two sources (see Braito et al. 2003a,b). 

At energies above 1 keV, the X-ray spectra appear to 
be dominated by a relatively flat power-law (PL) component 
in most of the sources. We tried to fit these overall 0.2- 
10 keV spectra with single-component models, either purely 
thermal or absorbed PL, as a zero-th order attempt. For all 
sources, these simple flts were rejected with high statistical 
significance. 

The XMiA-Newton spectral data are shown in Figs. 4 
and 5 and compared with model spectra. The bottom panel 
in each figures displays the ratio of the observed to predicted 
fiux as a function of energy. No emission lines are evident in 
these spectra, if we exclude Mkn 231 and IRAS19254-7245. 
The only appreciable feature is a fairly sharp, well charac- 
terized peak at E~ 0.7 — 1 keV, likely of thermal origin. 
We have modelled this component with a thermal emission 
template using the routine MEKAL in XSPEC (Mewe et al. 
1985), whose free parameters are the temperature kT and 
the normalization. The plasma metallicity has been set, for 
simplicity, to the solar value. This has no effect in our esti- 
mate of the the plasma temperature, the assumed metallicity 
essentially affecting only the normalization of the thermal 
spectrum. 
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Table 1. The source sample and summary of the XMM-JVetoton observations 



Name 


Coordinates 


z 




PN Net cxp. time 


Obs. Date 


Sequence 


Filter 








[iU L,Q\ 


n 1 

1 KSCC 1 


1 111011x11 / yCcir 1 






IRAS 12112+0305 


12 13 46.03 +02 48 41.5 


0.072 


3.8 


20.9 


2001 Dec 30 


0081340801 


Medium (PN) -Thin 


MKN231" 


12 56 14.16 +56 52 24.9 


0.042 


6.7 


19.8 


2001 Jun 06 


0081340201 


Medium 


IRAS 14348-1447 


14 37 38.30 -15 00 23.0 


0.082 


4.0 


17.5 


2002 Jul 29 


0081341401 


Medium 


IRAS 15250+3609 


15 26 59.43 +35 58 37.4 


0.055 


2.0 


17.8 


2002 Feb 22 


0081341101 


Thin 


IRAS 17208-0014 


17 23 22.03 -00 17 00.3 


0.043 


4.7 


15.7 


2002 Feb 19 


0081340601 


Medium 


IRAS 19254-7245*' 


19 31 21.46 -72 39 21.6 


0.062 


2.2 


18.4 


2001 Mar 21 


0081341001 


Thin 


IRAS 20100-4156 


20 13 29.75 -41 47 34.0 


0.129 


7.4 


18.1 


2001 Apr 21 


0081340501 


Medium 


IRAS 20551-4250 


20 58 27.05 -42 39 06.8 


0.043 


2.1 


16.1 


2001 Apr 21 


0081340401 


Thin 


IRAS 22491-1808 


22 51 49.26 -17 52 24.0 


0.078 


2.7 


21.7 


2001 May 24 


0081340901 


Medium 


IRAS 23128-5919 


23 15 47.00 -59 03 17.0 


0.044 


2.0 


9.4 


2002 Nov 19 


0081340301 


Medium 



Note: XMM-Newton observations reported in Braito et al. 2003b; 
XMM-Newton observations reported in Braito et al. 2003a 



In the rest of the present Section we discuss each source 
separately. The best-fit parameters are summarized in Table 
3. More detailed spectral decompositions are discussed in 
Sect. 5. 

4.1 IRAS 19254-7245, the Superantennae 

From the analysis of the ASCA data (Imanishi & Ucno 1999; 
Pappa, Georgantopoulos, & Stewart 2000), two alternative 
fits were proposed for this object: a) unabsorbed power law 
model with flat photon index (F ~ 1), and b) absorbed PL 
model with F ~ 1.8 and ^ 10^^ cm"^. No Fe-K lines 
were detected by ASCA. 

XMM-Newton has detected a strong X-ray flux from 
IRAS 19254-7245. As shown in Fig. 1, the emission is cen- 
tered in the southern nucleus of the interacting pair, while 
no hard X-ray emission is apparently related to the northern 
one. This is consistent with the results from observations at 
other wavelengths, which show that the southern nucleus, 
classified as a Seyfert 2, is always brighter than the nothern 
one (Melnick & Mirabel 1990). By taking into account the 
EPIC PSF, the X-ray source appears to be significantly ex- 
tended at low-energies (E< 2 keV). A detailed analysis of 
this complex source, resembling in many aspects the proto- 
typical type-II quasar NGC 6240, is reported in Braito et al. 
(2003a). 

A single unabsorbed power-law model is not a good 
description of the XMM-Newton data because it seriously 
misfits the spectrum at low energies. By adding a thermal 
component to their composite model, Braito et al. (2003a) 
have obtained a temperature kT=0.85lQ J keV for the ther- 
mal component and F = 1.84±0.20 (Nh = 4.7 x lO^^cm'^) 
for the hard PL component. The soft X-ray luminosity is 
~ 1.8 X 10"^ erg s~^ Braito et al. (2003a) have also found 
evidence for a Fe-K line at 6.4 keV, whose large equivalent 
width (EW~ 1.4 keV), together with the flat photon index, 
indicates that this object is a "Compton thick" AGN, in 
which the detected hard X-ray emission is due to a pure re- 
flected component plus a scattered one. Hence, the intrinsic 
AGN luminosity is much higher than the observed value of 
~ lO-*^ erg s-\ and likely > lO'** erg s'^ 



4.2 MKN 231 

Among local ULIRGs, Mkn231 is the most luminous object 
in the IR (Soifer et al. 1984) and one of the best studied 
at all wavelengths. Observed with many X-ray instruments 
(ROSAT, ASCA and Chandra), its X-ray properties still re- 
main puzzling. Evidences of a combined AGN and SB activ- 
ity emerged from the ROSAT and ASCA data (Imanishi and 
Ueno 1999; Turner 1999). However the flatness of the X-ray 
spectrum at energies above 2 keV and the lack of detection of 
any strong Fe lines (Maloney and Reynolds 2000) appeared 
unusual. These results have been confirmed with a recent 
Chandra observation (F = 1.3, FeK fine with EW<188 eV; 
see Gallagher et al. 2002). 

Our X}s.lM-Newton field containing the source is shown 
in Fig. 2: the source appears very bright and extended. This 
extension is significantly in excess of the PSF only for low- 
energy photons (E< 2 keV), with a total diameter of ~ 1 
arcminute (~ 50 kpc), while for higher energy photons the 
source spatial profile is consistent with the PSF. A detailed 
analysis of the XMM-A^euiton data on Mkn 231, plus those 
from a deep Beppo-SAX observation extending up to 80 keV, 
is reported in a companion paper by Braito et al. (2003b). 
Our results also confirm the very flat spectrum and the de- 
tection of a moderately intense (EW~200 eV) Fe-Ka line at 
6.4 keV. The luminosity of the hard component is ~ 3 x lO'*^ 
erg s~^. The presence of a highly obscured AGN is then con- 
firmed by the BeppoSAX PDS data, but our XMM-Newton 
observations suggest the concomitant presence of an impor- 
tant starburst component. 



4.3 IRAS 20551-4250 and IRAS 23128-5919 

IRAS 20551-4250 has been previously detected with ASCA 
(Misaki et al. 1999), but at a very low significance level. 
XMM-Newton detects it with a high enough S/N to allow us 
to perform a detailed spectral analysis. The source appears 
as unresolved by EPIC (see Figs. 2). 

A good fit (Fig. 5) is obtained only with a three com- 
ponents model: a thermal model (kT=0.68lg;Q7 keV), and 
a "leaky-absorber" continuum, including an absorbed plus 
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Table 2. Source Fluxes and Luminosities. 



JN SiTClQ 




'^0.5-2 




'-'2-10 




r(c) 

-^0.5-2 


r(c) 

-^2-10 




1'N+M()S1-A1()S2 


/■/■(/ rni 




ri-jj an " 


'] 


[<>'!, s 


[c,,y. -1] 


IRAS 12112+0305 


256+21 


1.2 X 10" 


'14 


1.5 X IQ- 


-14 


3.4 X IQ-*! 


3.6 X 10-*! 


MKN231 


2535 ± 51 


1.1 X 10- 


13 


6.5 X IQ- 


-13 


2.2 X 10-42 


6.4 X 1042 


IRAS 14348-1447 


297 + 24 


1.8 X IQ- 


-14 


2.0 X IQ- 


-14 


8.9 X 10"'! 


5.5 X 10*4 


IRAS 15250+3609 


433 ± 26 


2.1 X 10" 


-14 


2.4 X 10" 


-14 


3.0 X 10*1 


3.0 X 10*4 


IRAS 17208-0014 


401 ± 31 


2.3 X 10" 


-14 


4.0 X 10" 


-14 


6.0 X 10*1 


3.6 X 10*4 


IRAS 19254-7245 


926 ± 33 


5.3 X 10- 


-14 


2.3 X 10- 


-13 


1.8 X 10*2 


2.9 X 10*2 


IRAS 20100-4156 


101 ± 17 


5.4 X 10" 


-16 


1.9 X 10- 


-14 


1.2 X 10*2 


1.6 X 10*2 


IRAS 20551-4250 


942 ± 39 


5.5 X 10- 


-14 


1.4 X 10- 


-13 


4.2 X 10*2 


7.9 X 10*2 


IRAS 22491-1808 


232 ± 19 


1.0 X 10" 


-14 


6.4 X 10- 


-IB 


1.6 X 10*1 


1.7 X 10*4 


IRAS 23128-5919 


833 ± 34 


4.95 X 10- 


-14 


1.7 X 10- 


-13 


1.0 X 10*2 


1.8 X 10*2 



total net counts used in the spectral fits, 
observed fluxes. 

luminosities are corrected for galactic and intrinsic (Nh) absorption. 



a non-absorbed power law spectrum with the same pho- 
ton index. Wc found that a good fit can be obtained with 
r = l.Stl'l and Nh =8x W'^^cm-^. The 6.4 keV Fe-K line 
is not detected, but the upper limit on its equivalent width 
(< 1 keV) is consistent with the Nh value obtained from 
the fit. The intrinsic 2-10 keV luminosity of this object, cor- 
rected for absorption, is ~ 7.0 x 10*^ erg , while the soft 
X-ray luminosity is 4.2 x 10*^ erg s"^. The high X-ray 
luminosity and the properties of the X-ray spectrum clearly 
suggest the presence of an extinguished AGN. 

The XMM-Newton observations of IRAS 23128-5919 
have been strongly affected by high background (only 9.4 
ksec of good data are available), and this prevented us 
from attempting detailed spectral fits. The source appears 
as point-like and relatively bright. The spectrum is rea- 
sonably well reproduced by the sum of a thermal model 
and a simple power-law with no absorption. However, the 
power-law's spectral index in this solution is unphysically 
hard (F = 1.04 ± 0.13), suggesting the presence of an 
absorbed component. We then adopted the same "leaky- 
absorber" model as discussed for IRAS 20551-4250, and 
found best-fit parameters kT=0.65±0.15 keV, F = I.67I0 34 
and Nh = 6.9± 3.5 x 1022cm-2. The luminosities of the ther- 
mal and de-absorbed power-law components are ~ 1.5 x lO*'^ 
and ~ 2.7 x 10*2crg s^^. The large gas column density and 
power-law luminosity clearly indicate the presence of an 
AGN, whereas a spectroscopic IR study by Charmandaris 
et al. (2002) revealed the predominance of a starburst. Our 
subsequent analysis in Sect. 5.3 suggests that IRAS 23128- 
5919 is a transition object, with properties intermediate be- 
tween those of an AGN and a SB. 

4.4 IRAS 14348-1447, IRAS 15250-|-3609, and 
IRAS 17208-0014 

All these sources have intermediate-quality detections. 
While the latter two are consistent with a point source emis- 
sion at all energies, IRAS 14348-1447 appears to be resolved 
in a bow-like structure of total size -^30 arcsec in the North- 
South direction. The close-up image in Fig. 3 shows also the 



presence of another relatively bright (5'2-iofcev ~ 5 x 10~^® 
[erg/cm^ /s]) source located only 20 arcsec away in the SE 
direction, a position in which no optical counterpart is de- 
tectable in the POSS image. 

IRAS 14348-1447 and IRAS 17208-0014 show evi- 
dence of a single-temperature plasma emission with kT~ 
0.62±0.08 keV and kT~ 0.75±0.08 keV respectively. The 
spectrum of IRAS 15250+3609, instead, is unique among 
those from our spectral survey in requiring the presence of 
multi-temperature plasmas (a two-temperature solution in- 
dicates kT~ 0.26±0.09 keV and kT~ 0.64±0.15 keV). 

In addition to these thermal emissions we find a very 
significant hard X-ray component in all three sources. Ac- 
ceptable spectral fits have then been obtained by combin- 
ing a Mokal-modeled thermal emissions and an absorbed 
PL having F ~ 2.2, except for IRAS 15250+3609 whose 
photon- index is flatter (F ~ 1.2) (see Table 4). In all sources 
some absorption in excess of the galactic value is required 
(A'^H > lO^^cm-2). All three sources have X-ray luminosities 
in the range 10*^ - 10*2 erg s"^ 



4.5 IRAS 12112-1-0305, IRAS 20100-4156 and 
IRAS 22491-1808 

All these sources are detected at the lowest level of signifi- 
cance. This has prevented us from performing any detailed 
spectral analyses. We only tried on these sources two sepa- 
rate fits with a thermal model and an absorbed PL. For all 
three sources, a single thermal or a PL model are rejected 
by the data. Successful spectral fits require a combination 
of thermal and PL emissions (see Table 4). 

The intrinsic hard (2 — 10 keV) X-ray luminosities are 
1.7 X 10*^crg s-i for IRAS 22491-1808, ~ 1.6 x 10*2erg 
s"^ for IRAS 20100-4156, and 3.6 x 10*^erg s"^ for IRAS 
12112+0305. 
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Table 3. Best-Fit Spectral Parameters. 



Name 


model(=) 




kT 


thermal 


r 




T (d) 

Li'' ■ 








'kv\\ 


'■< '■// 








rn,. l; 


IRAS 12112+0305 


A 
B 


14.9/17 
16.3/18 


0.80«;- 


0.8 X 

1.1 X 


10^1 
10« 


1 92+"-33 
l-y^_0.55 

1.1* 


0.8+J-4 X 1021 

■/ 


6.3 X 104i 
6.1 X 1041 


MKN23l('') 




196 8/156 


n 07+0.09 


3.7 X 


10« 


1 30+°-^'' 


7 39+0 19 X 10^1 


8.3 X 1042 


TRAS 14S4S— 1447 


A 
B 


10.3/15 
13.7/16 


61+°-20 
0.64l°;r3 


3.5 X 
3.7 X 


10^1 

1041 


2 18+°-5i 
1.1* 


n 3+° '^ X 1022 
3 

0.5 j X 1021 


1.3 X 10*2 

1.1 X 1042 


IRAS 15250+3609(<') 


A 

■D 


13.9/18 

1/1 fl /I Q 

14. y/ lq 


64+°-^^ 

26+"-''** 
"•^"-0.07 

0-oo_o.io 
26+° °'^ 


2.3 X 
2.2 X 


10*1 
10« 


-I 9-1+0.62 
1-^1-0.66 

1 1 * 
l.i 


/ 
/ 


3.9 X 1041 


IRAS 17208-0014 


A 
B 


16.3/24 
17.1/24 


76 

"•'D -0.12 

74 +° " 


J-.D A 

1.3 X 


10*1 


2 26+1-1® 


1.1 +1-4 X 1022 

2.6 til X 10"' 


8.1 X 1041 

5.3 X 104i 


IRAS 19254-7245('') 


C 


48.1/56 


85+° " 


4.5 X 


1041 


1 84+°-®° 


0.47t°;g X 1022 


5.2 X 1042 


IRAS 20100-4156 


A 
B 


4.8/8 
4.7/8 


"■'^-0.32 
(, 7^+0.23 
"■'^-0.26 


4.1 X 
4.1 X 


10*1 

1041 


1.7* 

1.1* 


2.6+4-^ X 1022 
2.2^4:8 X 1022 


2.3 X 1042 

2.0 X 1042 


IRAS 20551-4250 


E(/) 


45.7/49 


66+°-^° 


3.0 X 


1041 


1 Q+0.20 


7.9t\l X 1023 


7.0 X 1042 


IRAS 22491-1808 


A 
B 


22.9/21 
26.7/23 


69+°-^^ 
u.oy_g 21 

0.661°;^^ 


1.5 X 
1.8 X 


10*1 

1041 


1 98+°-^^ 

1.1* 


/ 
/ 


3.6 X 1041 
3.8 X 104i 


IRAS 23128-5919 


E(/) 


32.7/44 




1.5 X 


1041 


1 67+°-°^ 
^•"'-0.34 


6.9t4-7 X 1022 


2.7 X 1042 



(a) : For this source an acceptable fit requires a two-temperature plasma. 
: XMM-JVewtoD observations reported in Braito et al. (2003a and 2003b) 
Model A: Thermal emission plus an absorbed power-law model; Model B: thermal emission plus eutoff power-law; Model C; 
Thermal emission plus a Compton reflected and scattered countinum (see Braito et al. 2003a); Model D: two thermal emission plus a 
"leaJfy-absorber" ; Model E: thermal emission plus a "leaky- absorber" continuum. 
Luminosities are between 0.5 and 10 keV. 
W): The covering factor for the "leaky-absorber" continuum are 9%, 5% and 54% for MKN231, IRAS 20551-4250 and IRAS 

23128-5919 respectively 
N.B. The symbol * indicate that the parameter has been kept fixed. 



5 THE SPECTRAL ANALYSIS 

The X-ray spectral properties of the ULIRG sample appear 
to be fairly uniform. The spectra are well reproduced by the 
combination of a soft thermal and a hard PL components. 
For a physical interpretation of the ULIRG phenomenon, 
the implications of these results are briefly discussed in this 
and the next Section. A more detailed comparison of these 
XMM-Newton spectra with model predictions is deferred to 
Persic et al. (2003a). 

5.1 Evidence for AGN-dominated emission 

One of the motivations for the present survey hcis been the 
opportunity offered by hard X-ray data to constrain the rel- 



ative AGN and staxburst contributions to the source ac- 
tivity. A dominant AGN contribution should be detectable 

in X-rays as cither a) an high- luminosity X-ray emission, 
L2-iofceV > 10''^ erg/s\ or h) highly extinguished hard X- 
ray components with Njf > 10^^ cm~^, as revealed by very 
flat or inverted hard X-ray spectra; or cj Fe-K complexes 
at ~6.4 keV with large Equivalent Widths {EW ~1 keV), 
corresponding to Iron fluorescent emission by cold molecular 
material illuminated by the AGN's energetic radiation field. 

We found such evidence for an AGN spectrum in 4 out 
of the 10 sources observed. The most obvious AGN is IRAS 
19254-7245, which has both a strong Fc-K line and a very flat 
(probably scattered) hard continuum, implying an intrinsic 
luminosity L2-10 > 10** erg/s (Braito et al. 2003a). 
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Figure 4. XM'M- Newton spectra in energy units and ratios of data to the best-fit model values as a function of energy, for the sources 
IRAS12112+0305, IRAS14348-1447, IRAS15250+3609, and IRAS17208-0014, from top left to bottom right. 



10 



The ULIRG Mkn 231 has been clearly identified from 
optical/IR observations as an AGN-dominated source, as in- 
ferred in particular from the flat, almost power-law, mid-IR 
continuum observed by Genzel et al. (1998) and the low 
EW of the PAH emissions. The source is also classified as 
a Broad Absorption Line quasar from optical spectroscopy 
(Smith et al.l995). The XMM-A''euitora observations reveal a 
moderate luminosity X-ray source with a very flat spectrum 
(r ~ 0.9) between 1 and 10 keV, while Beppo-SAX data 
reported by Braito et al. (2003b) reveal an highly extin- 
guished {Nh > 10^*cm~^) component above 10 keV. These 
data fully confirm the classification of Mkn 231 as a dusty 
and extinguished AGN. 

IRAS 20551-4250 is a third ULIRG with evidence from 
the XMM-A^euiton spectrum of harboring a high-column- 
density (Nh ~ 10^*cm~^) extinguished AGN, and a sim- 
ilar evidence for the presence of a low- luminosity AGN was 
found in IRAS 23128-5919. Neither these two sources nor 
IRAS 19254-7245 showed evidence for AGN activity from 
infrared spectroscopy (Genzel et al. 1998), which illustrates 
the useful complementarity between the IR and the X-ray 
approach. A detailed spectrophotometric study of Superan- 



tennae by Berta et al. (2003) has found evidence, however, 
for a combined AGN and starburst emissions with compa- 
rable IR bolometric luminosities. 

Further indications about the relative AGN/starburst 
contributions will come from the comparison of soft X-ray 
thermal and hard X-ray power-law emissions discussed in 
Sect. 5.3. 



5.2 Starburst emission 

The spectra of local SB galaxies in the 0.5-10 keV band can 
be described as a combination of warm thermal emission 
(with typically kT ~ 0.6 — 0.8 keV) dominating at energies 
< 1 keV, plus a harder spectrum producing the bulk of the 
2-10 keV fiux. 

The warm thermal component is interpreted as origi- 
nating from the boundary region between the hot, tenuous 
outgoing galactic wind and the cool, dense ISM (e.g., Strick- 
land et al. 2000). 

The harder component has various interpretations in 



© 0000 RAS, MNRAS 000, 1-?? 



10 A. Franceschini, V. Braito, M. Persic, et al. 



IRAS 20100-4156 



IRAS 20551-4250 



Unfolded Spectrum 



Unfolded Spectrum 



> I 

a a 



1 ' ' 

: ^ ' __; |_i ^ 

i — 1 — 1 — 1 — 1 — II 1 1 1 — i 


i 

~ *" 1 1 H 


o 
in 

>. '° 
m « 

(D 
1 

o 

a 
O 


1 ' ' 

1 ! 

1 1 i 1 1 


1 

1 1 : 







m 

o ^ 
o 





0.5 



Unfolded Spectrum 



1 2 
channel energy (keV) 

IRAS 22491-1808 



0.2 



0.5 1 2 

channel energy (keV) 

IRAS 19254-7245 



10 




* Unfolded Spectrum 
o 



> 



I 

o 

O CM 



— 1 1 1 1 1 — r 




1 2 
channel energy (keV) 

MKN 231 



7 Unfolded Spectrum 



0.5 



Unfolded Spectrum 



1 2 
channel energy (keV) 

IRAS 23128-5919 



10 



« O 



^ 



— 1 — 



1 .L-=^^^yi>i^ 



E 



5 

^ (O 




0.2 



0.5 



1 2 
channel energy (keV) 



10 



0.5 



10 



channel energy (keV) 



Figure 5. 'KMM-Newton spectra in physical units and ratios of data to the best-fit model, for the sources IRAS20100-4156, IRAS 
20551-4250, IRAS 22491-1808, IRAS 19254-7245, Mkn 231 and IRAS 23128-5919. 



the literature ^, either in terms of a very hot (kT>5 keV) 



^ Various authors have analysed X-ray data on specific nearby 
starbursts: Ptak et al. (1997) and Cappi et al. (1999) for NGC 253 



and M 82; Okada et al. (1997) for M83; Delia Ceca et al. (1996) 
for NGC 1569; Delia Ceca et al. (1999) for NGC 2146; Moran 
et al. (1999) for NGC 3256; Delia Ceca et al 1997 for NGC449; 



© 0000 RAS, MNRAS 000, 1-?? 



An XMM-Newton hard X-ray survey of ULIRGs 11 



„o.i 



0.01 



(a) 



2a491{}- 



.14348 



19254 

o 

23128 

o 



MKN231 

o 



10« 104; 
K.,J0.5 -10 keV) [erg s-'] 



10« 







"A 




/ 

/ 

/ 

/ 

/ 




- 


/ 

y 

// 20,006'''" 


- 




MKN23r 







14348 


/ ^ 


20551 


•E} / 







15250 / 




H} / 


23128 




22'191-t!- / 
^ I72O8 







/ 






7 / 12112 






/ -K- 






: / 






- / 






/ 




* - 



10*2 W' 

i„w (2-10 keV) [erg s-'] 



10*2 



(c) 



19254 O MKN2310 
-» 

20551 O {J.l'134a 
1 5250 

„ 17208 
23128 O ^ 



LriR [erg s-1] 



10" - 



ilO« 



10*1 



(d) 



MKN23 1 

o O 

20551 



O 

19254 



Q23138 



20100 



14348/ 

^■^^'^17208 

15250/ 
/ 

/ HJ 22491 



1048 
LpiE [erg s-i] 



J- 



(e) 



■a-15250 
■K-14348 
<J22491-»a0100 
{}17208 

{J-12112 



19254 

^ 20551 

O 

23128 
OqMKN231 



Figure 6. Luminosity-luminosity plots of various emission components identified in the XMM-Afeiuion spectra of the ULIRG sample. 
Starred symbols refer to sources that we classify as dominated by starbursts, open circles to AGN-dominated sources. Panel (a): luminosity 
ratio between the thermal and the power-law (2-10 keV) component after correction for intrinsic absorption, versus the bolometric X-ray 
luminosity measured in the total 0.5-10 keV band. Panel (h): luminosity-luminosity plot of the thermal versus PL components. Panels 
(c) and (d): luminosities of the thermal and PL components against the bolometric far-IR luminosities Li?/^. Panel {e): the ratios of the 
thermal and PL luminosities to the bolometric far-IR luminosity Li?/^. Dashed lines indicate proportionality relations. 
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Table 4. Observed ULIRGs Classification 



N^&rtie 


Optical 


Mid-IR 






\^) 




y^) 


i-tlAo Izliz+UoUo 


T 

ij 




oB 


Mkn 231 


AGN-1 


AGN 


AGN 


IRAS 14348-1447 


L 


SB 


SB 


IRAS 15250+3609 


L 


SB 


SB 


IRAS 17208-0014 


HII 


SB 


SB 


IRAS 19254^7245 


AGN-2 


SB 


AGN 


IRAS 20100-4156 


HII 


SB 


SB/AGN 


IRAS 20551-4250 


HII 


SB 


AGN 


IRAS 22491-1808 


HII 


SB 


SB 


IRAS 23128-5919 


HII 


SB 


SB/AGN 



Note. — Col.(l): Object name. Col.(2): Optical classification: 
Liner (L). AGN or HII region (Lutz et al. 1999; Veilleux et al. 
1999). Col.(S): Mid Infrared classification based on ISO spec- 
troscopy Starburst (SB) or AGN (Genzel et al. 1998). Col.(4): 
X-ray Classification. 

thermal or a F ~ 2 PL model, but no definitive conclusions 
were reached, as the thermal and PL fits are, in general, 
similarly successful. 

Recently, Persic & Rephaeli (2002; see also David, Jones 
& Forman 1992) have developed a detailed quantitative 
model of synthetic X-ray spectra of SB galaxies, based on 
evolutionary populations of galactic stars and adopting tem- 
plate X-ray spectra for the relevant emission processes. They 
suggested that high- and low-mass X-ray binaries (HMXBs, 
LMXBs) contribute most of the 2-15 keV spectrum, in the 
absence of AGN emission. Both categories of X-ray binaries 
have spectra that can be described as variously cutoff PLs 
(White et al. 1983; Christian & Swank 1997). For example, 
in the case of one isolated episode of intense star formation, 
the 2-10 keV emission of the starburst would be dominated 
by HMXBs, and would be described as a F = 1.2 power 
law. In more moderate starbursts (such as those observed 
in local SB galaxies), in addition to the HMXBs related to 
the SB proper, also the LMXB population in the underlying 
disk is expected to provide important contributions to the 
X-ray emission. In the case of a mix of HMXBs and LMXBs 
of various luminosities with Galactic proportions, the 2-10 
keV emission can be described as a cutofT PL of the form 

exp~^^'°^ , with photon index F ~ 1.1 and cutoff en- 
ergy feT ~ 8 keV (Persic & Rephaeh 2002). This spectrum 
is to be corrected for the intrinsic absorption with variable 
column density Njf. 

5.3 Spectral decompositions 

Following the above guidelines, we have tried to fit the spec- 
tra of our ULIRGs with various astrophysically motivated 
composites. The first clear feature to consider is the ubiq- 
uitous low-energy excess emission between 0.5 and 1 keV. 
One possibility to explain it would be a reflected component 
by photoionized g found for example in the Seyfert 2 

NGC 1068 (Kinkhabwala et al. 2002, Brinkman et al. 2002). 

Zezas et al. (1998) and Delia Ceca et al. (2002) for NGC 3310, 
NGC 3690 (Arp 299). See also Dahlem et al. (1998) for a review. 



We have tried to reproduce the source spectra with such 
an emission, but in all cases have failed to obtain acceptable 
fits. In particular, for IRAS 19254-7245 the Fe-K Une at 6.49 
keV detected by Braito et al. (2003a) implies the presence 
of an essentially cold neutral gas, which is unable to explain 
the low-energy emission features. For MKN 231, Braito et al. 
(2003b) attempt to explain it entirely with an ionized reflec- 
tion, which provides however a very poor overall fit, while 
a thermal ionized gas emission is required by the data. As 
for IRAS 20551-4250, an attempt to fit the low-energy spec- 
trum with photo-ionized emission is similarly unsuccessful, 
as it would require an unphysically large value of the X-ray 
spectral index F and would generate a quite poor overall fit. 

Also in all other sample's sources the photo-ionized 
AGN emission model would leave significant residuals in all 
cases, which implies that coUisionally ionized hot thermal 
gas is in any case required. We have modelled this ubiqui- 
tous low-energy thermal component with the Mekal/XSPEC 
code at constant (solar) metallicity (for the source IRAS 
15250-1-3609 a two- temperature component was required). 

Then we have added a PL with photon index F, pho- 
toelectrically absorbed through a column density Nh to 
describe the hard X-ray excess that shows up at energies 
> 2 keV. For IRAS 20551-4250 we have used the "leaky ab- 
sorber" model (Sect. 4.3), while for Mkn 231 we have not 
considered in the present analysis the highly-extinguished 
(Nh > lO'^^cm,^^) component found in Beppo-SAX data by 
Braito et al. (2003b). Alternatively, for all the sources we 
have attempted also to fit the high-energy data with an X- 
ray binary model with F = 1.1, cutoff energy kT ~ 8 keV 
and variable low-energy absorbing column density Njy (see 
Sect. 5.2). 

The results of these spectral decompositions are re- 
ported in Table 4 and shown in Figs. 4 and 5. In Fig. 6 
we explore possible relationships between the best-fit pa- 
rameters. Fig.6(aj is a plot of the ratio between the lumi- 
nosities of the thermal and PL components (after correction 
for intrinsic absorption) , as a function of the X-ray luminos- 
ity measured in the total 0.5-10 keV band. SB-dominated 
ULIRGs (staxred symbols) occupy a fairly well defined 
region with Ltotai ~ W^^erg/s and 

IRAS 12112+0305 appears to have a very low value 
of Lthermai Compared with other starbursts. The AGN- 
dominated sources have enhanced X-ray luminosities with 
^thermal /^PL < 0.1: as expected, these AGN ULIRGs have 
an energetically-dominant power-law component. In the fol- 
lowing we use the symbol Lpi to indicate the luminosity 
derived from the fit with a generic PL, rather than the one 
corresponding to the best-fit X-ray binary model (in any 
case the difference between the two is a few tens percent at 
most, see Table 4). 

Fig. Q(b) compares the luminosities of the thermal and 
PL components: a fairly clear correlation is apparent, with 
the AGN population occupying the high-luminosity end in 
both quantities. Again the AGN-dominatcd sources show an 
excess of the PL component luminosity compared with that 
of the thermal plasma. 

Fig. 6(c) plots the luminosities of the thermal com- 
ponent Ifthermai agaiust the bolometric far-IR luminosities 
Lf/a, while Fig. G(d) does the same for the luminosities of 
the PL component. If we consider the starburst-dominated 
population, the Lpl and Lfir luminosities appear both rear- 
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sonably correlated, while the AGN-dominatcd sources show 
similar Lfir values but largely enhanced X-ray PL emission. 
On the contrary, we find a remarkable lack of correlation be- 
tween Lfir and the emissivity of the thermal component. 

Finally, panel (3(e) shows that the SB-dominated sources 
have fairly well defined Lx I^fir ratios for both the thermal 
and PL components. The AGN sources display the usual 
excess of PL emission. On this regard, the sources IRAS 
20100-4156 and 23128-5919 tend to occupy positions in the 
various plots which are intermediate between SB and AGN 
dominated populations, and may have genuinely intermedi- 
ate properties. 



6 DISCUSSION 

Altogether, in spite of the limited statistics of the sample, 
there are some remarkable regularities emerging from the 
spectral analysis of the XMM-ATewton data. 

The first one is that hot thermal plasma emission, likely 
a SB signature, is present in all objects of the sample. We 
find, however, a remarkable lack of correlation between ther- 
mal emission and far-IR luminosity (which is supposed to be 
a good indicator of the rate of star formation in the galaxy). 
Our conclusion is that, although this plasma emission should 
be traced back to a galactic wind triggered by young and ex- 
ploding stars, the process may be ruled by various additional 
parameters, like the density and pressure of the surrounding 
ISM. According to the study by Strickland & Stevens (2000), 
the soft X-ray emission from SB-driven galactic winds comes 
from a low filling-factor (< 2%) gas, which contains only a 
small fraction (< 10%) of the mass and energy of the wind. 
Soft X-ray observations therefore do not probe the gas that 
contains the bulk of the energy, mass or metals in the out- 
flow, while the bulk of the hot plasma radiates and cools on 
timescales longer than that of the starburst evolution. 

To test the possible effect of the long cooling time, we 
have attempted to compare the thermal X-ray emissivity 
with the starburst ages. To this end, since all optical coun- 
terparts of our sources show double or multiple nuclei (sup- 
posed in a merging system), we have adopted as a rough 
estimate of the merger stage the separation of the nuclear 
components, and investigated possible relations between the 
ratio of the thermal X-ray to IR luminosities against nuclear 
separation. We found essentially no correlations in these 
plots, which suggests that various parameters in addition 
to the star-formation rate (SFR), like the pressure of the 
surrounding medium, should affect the thermal plasma emis- 
sion. In any case, the soft thermal emission is a quite poor 
tracer of the ongoing SFR. 

On the contrary, the PL emission fitting the hard X- 
ray excesses in the spectra of our SB-dominated ULIRGs 
seems well correlated with Lf/«. We understand this as an 
effect of both quantities being closely linked to the newly 
formed stellar populations, in one case (Lf/h) due to dust 
reprocessing of the UV flux by young massive stars, in the 
other (Lp_l) to the number of HMXBs which are a subset of 
the young stellar population (see below) . 



6.1 Stellar contribution to the X-ray emission 

We now check the hypothesis that, in ULIRGs not domi- 
nated by an AGN, the hard X-ray 2-10 keV emission may 
be mainly due to luminous (I/2-10 > lO^'^ erg s~^) HMXBs. 
To this aim, we compare the rate of star-formation (SFR) 
estimated from the X-ray flux with that traced by the FIR, 
under the assumption that both the X-ray emitting HMXBs 
and the OB stars heating the dust responsible for the FIR 
flux are short-lived, in which case they can both be used as 
indicators of the ongoing SFR. Our definition of the SFR 
assumes a standard Salpeter stellar IMF between 0.1 and 
100 Mq. 

We proceed £is follows. Assuming a mean HMXB lu- 
minosity of L2-10 = 5 X 10^'' erg s~^ (e.g., White et al. 
1983), we first estimate the number of HMXBs from the 2- 
10 keV luminosity for each sample object. We consider that 
our Galaxy, hosting ~ 50 bright HMXBs (Iben et al. 1995, 
and references therein), has a SFR of ~ ZMq yr~^ (e-g-, 
Matteucci 2002). From the number of HMXBs we then es- 
timate the corresponding SFRx-ray 

Then from the FIR luminosity, we calculate (Kennicutt 
1998) 

SFReY^ ~ LFm/(2.2 x lO^^erff s-')Moj/r-\ 

We plot the two independent estimates of the SFR in Fig. 
7: this shows that the two are in quite good agreement for 
these non-AGN-dominated objects (filled circles), whereas 
the SFRx -ray valucs are clearly in excess for the AGN- 
dominated sources (empty squares). Given the linear rela- 
tion between SFR, far-IR and X-ray luminosities, this is ob- 
viously nothing else than the Lpir — L2-iokcV plot of Fig. 
6d, recast in different units. 

We have further checked this X-ray to FIR relationship 
on a sample of well-known moderate-to-low luminosity local 
SB galaxies, shown in Fig. 7 as crosses (see figure's caption 
for references). While comparing these with ULIRGs, and 
following our previous discussion in Sect. 5.2, we bear in 
mind one relevant difference between the two classes: the 
star-forming activity in ULIRGs is very intense, short-lived, 
and dominates the whole galaxy, so the hard X-ray flux is 
produced mainly by HMXBs. In lower SFR starburst galax- 
ies the X-ray flux is produced by a mix of HMXBs and 
LMXBs, coming respectively from SB regions and from the 
underlying quiescent disk, in comparable proportions. In 
Fig. 7 we have quantified this effect by the parameter /, 
the fraction of 2-10 keV emission attributable to HMXBs, 
setting / = 1 for ULIRGs and / ~ 0.25 for SB's (Persic & 
Rephaefi 2002). 

These results confirm that, for ULIRGs without domi- 
nant AGN components, the 2-10 keV flux from ULIRGs is 
a good SFR indicator: 



SFR: 



'X — ray 



(io55i;7o 



(1) 



for L2-i() > W^^erg s^^ . For lower luminosity SB galaxies 
(Z/2-10 < W^^erg s~^), the relation of the X-ray flux to 
the SFR is less straightforward, due to the contribution of 
long-lived LMXBs. In this case the relation becomes 

" - (o4) (4xl£C.-^) ^^^'^"'- 
These results are in fair agreement with those found by 
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Figure 7. - The SFR estimated from X-ray emission versus that estimated from FIR emission. Empty squares and filled dots represent, 
respectively, AGN-dominated and non- AGN-dominated ULIRGs (from our sample) , while crosses denote SB galaxies (from the literature) . 

We used a representative SB sample (sec section 5.2), for which: a) distances and FIR luminosities come from Shapley et al. 2001 (M82, 
M83, NGC 253, NGC 1569, NGC 3079, NGC 3310, NGC 3628, NGC 3690, NGC 4631), Dahlcm ct al. 1998 (NGC 55, NGC 4666), Delia 
Ceca et al. 1999 (NGC 2146), and Moran ct al. 1999 (NGC 3256); and b) 2-10 keV fluxes come from Dahlcm et al. 1998 (M82, NGC 
55, NGC 253, NGC 3079, NGC 3628, NGC 4631), Okada et al. 1997 (M83), Delia Ceca et al. 1996 (NGC 1569), Delia Ceca ct al. 1999 
(NGC 2146), Moran et al. 1999 (NGC 3256), Zezas et al. 1998 (NGC 3310, NGC 3690), and Persic et al. 2003b (NGC 4666). The dotted 
hue marks the relation SFRx-ray = SFRpiR. The source with SFR^ ~750 and SFKfir ~50MQ/yr corresponds to NGC 3690 (ARP 
299), which was found by Delia Ceca et al. (2002) to contain an absorbed AGN. 



Ranalli ct al. (2002) and Gilfanov ct al. (2003), if we con- 
sider the different definitions of the SFR. Note that the only 
SB galaxy lying far out the one-to-one relation in Fig. 7 is 
the crossed symbol with SFRx-ray ~ 750 Al^yr'^ corre- 
sponding to the galaxy Arp 299, which was recently proven 
to host a luminous obscured AGN (Delia Ceca et al. 2002). 

Finally, our XMM-Newton data below ~ 0.5 keV have 
in all cases a very low S/N, which does not allow us to con- 
strain the origin of this part of the ULIRG spectrum. In some 
instances (IRAS12112-F0305 and IRAS 22491-1808) our for- 
mal fit implies a dominance of X-ray binary emission at such 



low energies, which may not be physical due to the fact that 
binary spectra typically show photoelectric absorption. To 
address this problem with more elaborated model spectra 
would need higher spectral resolution and S/N data. 

6.2 The AGN contribution to the ULIRG 
phenomenon 

Altogether, we find in 3 of the 10 sample sources (Mkn 231, 
IRAS 19254-7245 and IRAS 20551-4250) various indepen- 
dent evidences for the presence of an absorbed AGN dom- 



© 0000 RAS, MNRAS 000, 1-?? 



An XMM-Newton hard X-ray survey of ULIRGs 15 



inating the XMM-A'^euiton X-ray spectrum. They show, in 
particular, values of the 2-10 keV luminosity (after absorp- 
tion correction) quite in excess compared with the thermal 
plasma luminosities. This evidence is confirmed in all three 
objects by X-ray spectral features like a strong Fe-K 6.4 
kcV lino and a very flat or inverted hard X-ray spectrum in- 
dicative of high-column density {Nh > 10^^ cm^^) circum- 
nuclear material. 

The source IRAS 23128-5919 displays a hard spectrum 
and large X-ray luminosity which suggest an AGN contribu- 
tion. Both IRAS 23128-5919 and IRAS 20100-4156 display 
intermediate values of the parameters in Fig. 6 between the 
AGN and SB-dominated sources, which might indicate that 
an AGN contribution may be present. 

It should be noticed that, among X-ray AGNs, IRAS 
19254-7245, IRAS 20551-4250 and IRAS 23128-5919 did not 
show any AGN signatures from spectroscopy of IR coronal 
lines (Genzel et al. 1998). Consider however that, as shown 
in Fig. 6e, the X-ray luminosity is only a tiny fraction (~ 
0.01% and 0.1% for SBs and AGNs, respectively) of the 
bolometric one. Then even for bona-Gde X-ray AGNs, the 
IR spectrum may well be dominated by the starburst, hence 
explaining the results of the IR spectroscopy. Although the 
lack of AGN signatures in the IR spectra of those sources 
might alternatively be explained as dust obscuration up to 
30 fim, the evidence for staxburst components in the X-ray 
spectra suggests that indeed an important fraction of the 
bolometric IR flux is likely of stellar origin (see previous 
Sect, and Fig. 7). 

Though limited by the small statistics of our ULIRG 
sample, our conclusion is that for a majority of these sources 
the high rate of star formation indicated by the large far-IR 
luminosity can account for most of the X-ray emission. The 
X-ray spectra of roughly half of the ULIRGs show evidence 
of AGN contributions on top of an, always present, starburst 
component. 



7 CONCLUSION 

We have devoted a large observing program with XMM- 
Newton to survey the hard X-ray properties of a complete 
and representative sample of Ultra-Luminous IR Galaxies, 
as a way to probe deeply into the heavily extinguished cores 
of this still physically elusive class of sources. 

All the 10 observed ULIRGs have been detected with 
high statistical significance by XMM-Newton, but their 2-10 
keV fiuxes turned out to be rather faint on average. 

The X-ray emission appears to be extended on a scale 
of -30 kpc for Mkn 231 and IRAS 19254-7245, possibly 
evidence of galactic superwinds. In these same sources, in 
IRAS 20551-4250 and IRAS 23128-5919 we find evidence for 
the presence of hidden AGNs, while a minor AGN contribu- 
tion may be suspected also in IRAS 20100-4156. A strong 
Fc-K line (EW~ 2 keV) in the X-ray spectrum of IRAS 
19254-7245 and a much weaker line in Mkn 231 (EW~ 0.2 
keV) are also detected, suggestive of deeply buried type-II 
quasars (Braito et al. 2003a,b). For the other sources, hence 
for roughly half of the ULIRG sample, the X-ray luminosi- 
ties and spectral shapes are mostly consistent with hot ther- 
mal plcisma and X-ray binary emissions of mainly starburst 
origin. 



We have analysed the X.MM-Newton spectra in terms of 
various physical components, that is thermal plasma emis- 
sion and hard X-ray power-laws due to either X-ray binaries 
or true AGNs. We have found interesting regularities in the 
source parameters. Thermal plasma emission, the signature 
of a starburst component and dominating the spectra be- 
tween 0.5 and 1 keV, is present in all the sample objects. 
The alternative interpretation of this soft X-ray component 
as due to a reflected emission by AGN-photoionized gas is 
mostly excluded by our analysis. 

This thermal emission is quite unrelated with the far- 
IR luminosity (a tracer of the ongoing star formation rate) , 
from which we infer that other parameters should determine 
it in addition to the rate of star formation. 

On the contrary, the X-ray binary power-law emis- 
sion fitting the hard X-ray component in the starburst- 
dominated ULIRGs is correlated with Lfir, both quanti- 
ties tracing the number of young stars in the galaxy and 
measuring the ongoing SFR. 

We fully confirm with these data the composite nature 
of ULIRGs as a class, with indications for a predominance 
of the starburst over the AGN phenomenon even when ob- 
served in hard X-rays. 
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